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We propose a method for controllable generation of non-local entangled pairs using spinor atoms 
loaded in an optical superlattice. Our scheme iteratively increases the distance between entangled 
atoms by controlling the coupling between the double wells. When implemented in a finite linear 
chain of 2N atoms, it creates a triplet valence bond state with large persistency of entanglement (of 
the order of N) . We also study the non-equilibrium dynamics of the one-dimensional ferromagnetic 
Heisenberg Hamiltonian and show that the time evolution of a state of decoupled triplets on each 
double well leads to the formation of a highly entangled state where short-distance antiferromagnetic 
correlations coexist with longer-distance ferromagnetic ones. We present methods for detection and 
characterization of the various dynamically generated states. These ideas are a step forward towards 
the use of atoms trapped by light as quantum information processors and quantum simulators. 



I. INTRODUCTION 

The generation and manipulation of entanglement 
have been identified as important requirements for quan- 
tum teleportation [1], quantum information processing 
[2] and quantum communication [3]. Engineering long- 
ranged entangled pairs in optical lattices can also have 
fundamental implications in the context of quantum mag- 
netism. For example, many frustrated spin states such as 
spin liquid states correspond to coherent superpositions 
of spin singlet states [4]. 

Recent experiments have made progress towards gen- 
erating multiparticlc entanglement among an ensemble of 
atoms confined in optical lattices by using controlled col- 
lisions between individual neighboring atoms [5]. How- 
ever, the generation of long-distance pair entanglement 
in systems with short-range interaction between particles 
(such as optical lattices) is not an easy task. In recent 
proposals long-distance EPR pairs [1] are generated by 
first creating an entangled pair of quantum particles in 
one location and then physically transporting one mem- 
ber of the pair to another location [6]. However, decoher- 
ence during the transport reduces the quality (fidelity) of 
the entanglement. 

Our approach is based on coherent manipulations of 
triplet or singlet pairs of ultra-cold atoms loaded in an ar- 
ray of double-well potentials called superlattice [7, 8, 9]. 
These manipulations, applied to isolated double wells, 
were used for the recent observation of supcrcxchangc 
interactions in optical lattices [10, 11, 12]. Here we gen- 
eralize these approaches to study the many-body dynam- 
ics that arises when coupling between the double wells 
is allowed for. We propose various schemes that result 
in controllable generation of multiparticle entanglement. 
Specifically, we first discuss a protocol that creates from 
a system of spinor bosonic atoms initially prepared as 
an array of triplet (singlet) pairs on neighboring sites, 
an array of long-distance triplet (singlet) pairs across the 
lattice. The method consists of a simple iterative swap- 
ping procedure, performed by controlling the double- well 



barrier height (see Fig. 1), which enables parallel gener- 
ation of long-distance EPR pairs. 

We find that by combining the iterative swapping pro- 
cedure with the boundary effects always present in a fi- 
nite linear chain, one can engineer a state in which each 
atom located in the right half of the superlattice is en- 
tangled with an atom in the left half. This bipartition 
of the system into its left and right parts exhibits maxi- 
mal entanglement entropy. Similar procedures have been 
proposed for coherently transporting quantum informa- 
tion [13] and for creating bosonic cooper-like pairs [14] 
in optical lattices. Additionally, we show that the par- 
allel generation of an array of EPR pairs can be useful 
for efficient implementation of entanglement purification 
schemes [15], which aim to distill the few high-fidelity 
entangled pairs from the numerous low-fidelity ones. 

The swapping procedure described above is imple- 
mented in an array of decoupled double wells. An in- 
teresting question that naturally arises is what happens 
with the state if the double wells are no longer completely 
decoupled, but instead there exists a finite tunneling be- 
tween them. The resulting dynamics goes beyond the 
simple two-particle physics behind the swapping proce- 
dure and the experiments which control superexchange 
interactions [10]. The emerging state is the consequence 
of many-body dynamics of a global interacting Hamilto- 
nian and does not require manipulations on individually 
accessed atoms. This is a promising approach for creat- 
ing new magnetic phases without explicitly processing a 
quantum-computer protocol. Although we believe that 
the phenomena we discuss here are very general, to be 
specific we consider in this paper a one-dimensional chain 
and focus on the coherent evolution of the product state 
of triplets or singlets in each individual double well (Fig. 
la). These are dimerized states which break translational 
symmetry. This choice of initial states is motivated by 
the fact they can be prepared in experiments [10]. 

Our analysis shows that the time evolution of the 
triplet product state leads to the formation of a magnetic 
state with mixed correlations and a high degree of mul- 
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tiparticle entanglement, where short-range antiferromag- 
netic and long-range ferromagnetic correlations coexist. 
This state can be experimentally probed by measuring 
the singlet-triplet populations [11] and density-density 
correlations after time of flight [16]. We also find total 
(partial) restoration of the translational (rotational) sym- 
metry, which suggests that our final state has some type 
of spin liquid character. By this we mean a state with 
strong intrinsic fluctuations but no broken symmetries 
[4], what may be different from other definitions which 
are based on the topological order of the quantum state 
[17]. 

The time evolution of the initial singlet state also leads 
to the restoration of the translational symmetry and high 
multiparticlc entanglement but in this case we do not 
observe the strongly mixed correlations. The dynamic 
state has purely antiferromagnetic character, although 
with an unusual behavior of long-range correlations. 

The paper is organized as follows: After introducing 
in Sec. II the formalism and numerical techniques we use 
for our analysis, in Sec. Ill we describe the basic Hamil- 
tonian and its possible implementation in the context of 
recent experiments using optical superlattices. In Sec. IV 
we present the swapping procedure which we refer to as 
a single switch dynamics and in Sec. V we discuss the 
idea of iterative repetition of the switch as a means to 
generate long-distance entangled pairs. We also study 
possible ways to experimentally detect such long-range 
correlations. In Sec. VI we relax the isolated double- we 11 
constraint and allow for a finite coupling between the 
double wells. Specifically, we concentrate our analysis on 
the many-body dynamics that emerges when both the 
intra- and inter-well couplings are equal and study the 
coherent dynamics starting from both an initially pre- 
pared triplet product state and an initially prepared sin- 
glet product state. Finally, we present our conclusions in 
Sec. VII. 



II. THE FORMALISM 

The focus of this paper is twofold. On one hand we 
study experimentally relevant observables which can be 
used to detect and characterize the dynamics of cold 
atoms. On the other hand we analyze properties of en- 
tanglement in the system. The propagation and redistri- 
bution of entanglement are not only important from the 
quantum- information perspective, but can also help to 
understand the quasiparticle dynamics as demonstrated 
recently [18]. Such properties are best discussed in terms 
of the entanglement entropy which corresponds to the 
von Neumann entropy of the reduced density matrix with 
respect to a bipartition into two subsystems [19]. The 
entanglement entropy is defined as S = — tr(plog 2 (p)), 
where the reduced density matrix p — tr|^>(t))(^>(t)| is 
the trace over the states of either of the two subsys- 
tems. For the one-dimensional systems with open bound- 
ary conditions, we will study the entanglement entropy 



Si of a block of size I located at the edge of the chain. 
In the case of an infinite system we define S^ en (S%~ ) 
as the entropy of subsystems formed by partitioning the 
chain at an even (odd) bond. While any product state (a 
state that can be represented as a tensor product of two 
pure subsystem states) has zero entanglement entropy, 
maximally entanglement states at half bipartition have 
entanglement entropy of S = N. 

We use both numerical and analytic techniques to 
study the quantum dynamics. For the numerical treat- 
ment we adopt the time-evolving block decimation algo- 
rithm (TEBD) for finite [20, 21] and periodic infinite sys- 
tems [22], which uses a matrix-product state representa- 
tion and a Suzuki- Trotter decomposition of the evolution 
operator. It retains only states with the lowest weights in 
the reduced density matrix, keeping the number of states 
X (the dimension of the matrices) finite. Consequently, 
the wave- function of weakly entangled states can be han- 
dled efficiently, with the computation times of the order 
ofO(x 3 iV). 

During the time evolution \ has to be increased in or- 
der to reproduce the growing entanglement in the system. 
The accuracy of the method is estimated by varying both 
X and the Suzuki- Trotter slicing [23]. For short and in- 
termediate times the TEBD algorithm allows us to get 
very precise results, but at the moment when the en- 
tanglement entropy exceeds log2(x)i the matrix-product 
representation becomes no longer accurate. To deal with 
the evolution over long periods of time (t oo), we use 
exact diagonalization [24] techniques. Even though these 
techniques can only deal with systems with small number 
of lattice sites (up to 24 sites) and suffer from recurrence 
effects, they are relevant for realistic setups in ID exper- 
iments [25]. 



III. SETUP AND PROCEDURES 
A. Effective Hamiltonian 

We consider a system of 2N ultracold bosonic atoms 
with two relevant hypcrfine states, which we denote as 
t and I, confined within a double- well superlattice with 
the filling factor of 1. The latter can be experimentally 
implemented by superimposing two independent lattices 
one with twice the period of the other [9, 10]. 

In the deep barrier regime, the vibrational energy of 
each well, Jtwq, is the largest energy scale in the sys- 
tem and one can restrict the dynamics to the lowest vi- 
brational states. When restricted to the lowest band, 
there are three relevant energy scales: the intra-well hop- 
ping amplitude tj„, the inter- well hopping amplitude t out 
and the on-site interaction energy U . In the limit of 
large U » t in ,t out we are focusing on, the system is 
in the Mott insulating regime and the only populated 
states are the singly occupied ones. The spin dynam- 
ics is described by the following effective Hamiltonian, 
which takes into account the coupling between the differ- 



ent singly occupied states by virtual particle-hole excita- 
tions [10, 26, 27], 

H e ^ = -Ji ^ S 2 j ■ S 2j) +i - J 2 ^ S 2 j+i • Saj+2 , (1) 

3 3 

with J 1 = 4:tf n /U and J 2 = 4tg Ut /U. Since experimen- 
tally ti n and t otl t can be controlled independently [11] by 
adjusting the intensities of the laser beams that generate 
the superlattice, we will assume that both J\ and </ 2 can, 
in general, be time- dependent functions J\{i) 1 J 2 {t). Ad- 
ditionally, we note that even though for bosons the sign 
of the coupling constants is normally positive (ferromag- 
netic interactions), experimentally it is also possible to 
change the sign to be negative [10]. 



B. Initial state 

The starting point of our analysis is a system initially 
prepared in an array of triplet pairs on the neighboring 
sites of a double- well superlattice, 

w = o)> = n i*wi> . ( 2 ) 
3 

l*L-+i) = ^(lt)il4) J -+i + U),-|t)i+i). (3) 

This state has been recently realized in the laboratory 
[10]. In this experiment, after first preparing a Mott 
insulator with two bosonic atoms per double well, the 
atoms were transferred into a triplet state configuration 
by using spin-changing collisions [28]. 

For the following, it is convenient to characterize the 
initial state as a triplet valence bond state of length 1. 
Although this state is a ground state of the system of 
independent wells, it is not an eigenstate of a system of 
coupled wells. Therefore, changing the couplings Ji,2 at 
t > leads to a complicated correlated dynamics. The 
specific time evolution depends significantly on the ratio 
of the couplings J\ and J 2 . 



C. Switching procedures 

We consider and characterize in details three specific 
cases: 

1. Single switch: J x [t > 0] = 0, J 2 [t > 0] = J. 

2. Periodic switch: Ji[(2n + l)t s > t > 2nt s ] = 
0,J 2 [(2n + l)t s > t > 2nt s ] = J while Ji[(2n + 
2)t s > t > (2n + l)t s ] = J, J 2 [(2n + 2)t s > t > 
(2n + l)t s ] = with n = 0,1,2,... and switching 
time t s specified below. 

3. Homogeneous switch: J\\t > 0] = J 2 [t > 0] = J. 
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FIG. 1: a) The initial state in the superlattice corresponds 
to a product of triplets at adjacent sites, a is the lattice spac- 
ing, b) At time t = + the intra-well tunneling is suppressed 
and the inter- well tunneling is allowed, c) At t — t s the en- 
tanglement between adjacent pairs is redistributed between 
pairs of length 3. d) If the switching procedure is repeated, 
the entanglement propagates to atoms separated by 5 wells 
and after n switches by 2n + 1 wells. 



The Hamiltonian in the first two cases consists of decou- 
pled double wells and allows a simple analytical treat- 
ment (Sections IV ,V). The homogeneous switch involves 
the complicated many-body dynamics of the Heisenberg 
chain and will be analyzed using numerical tools (Sec- 
tion VI). 

It is convenient to introduce the bond operators [29] 
which create singlet and triplet pairs at different bonds: 

4,+il°> = l*ij+i) = ^(IT)iU>j+i-U>ilt>i+i) 1 
t£- +1 |o) = (4) 
= l*L-+i) = 7|(IT)ilT) J - +1 -U) J -UW), 

= l*5W = Tfd t>il T)i+i + 1 ih+i) 

(|0) denotes the state with no atoms). These operators 
satisfy bosonic commutation relations and the constraint 

E ^+1^+1+4+1^ = 1 > ( 5 ) 



4 



which follows from the completeness of the Hilbert space 
of states of an individual double well. We start our anal- 
ysis by studying the single switch dynamics. 



IV. SINGLE SWITCH 

In the case J\ = and J 2 = J, the evolution operator 
U(t) = e~ ltE -l n can be written analytically. It is given by 



U odd (t) = 



iAt 



*sin(^)P 2 j+i,2i+2] , (6) 



where P2j+i,2j+2 = (1/2) + 2S 2j+ i ■ S 2j +2 is the swap 
operator (interchanges the spins) at sites 2j + 1, 2j + 2 
and A is an irrelevant phase factor equal to —J(N + 
2) /(Ah). From Eq. (6) it is clear that the evolution is 
periodic with the period 




FIG. 2: Entanglement entropy for the single switch and for 
the periodic switch (numerical result, TEBD), t s = irh/J. We 
used an infinite lattice for the former and one with 2N = 10 
for the latter and calculated the entanglement entropy for half 
of the chain. While for the single switch the period is 2t 3 , for 
the periodic switch the initial state is recovered after t = 
2Nt s . The single switch is well described by expression (9). 



T = 2t s , 
Jt s = irh. (7) 

At times t = (2n+l)t s , the evolution operator reduces 
to a product of the swap operators which, upon acting on 
the initial state, distribute the entanglement from atoms 
at sites (2j, 2 j + 1) to atoms at (2j + 1, 2j + 4), leading 
to the formation of a quantum state with valence bond 
length equal to 3 (see Fig.l). 

i^=*.> (i) =ni^+iw- ( g ) 

The effect of this redistribution on the entanglement 
entropy is shown on Fig. 2. We observe that while for odd 
bipartitions the entanglement entropy oscillates between 
and 2, for even bipartitions S^ en remains constantly 1. 
This is consistent with the fact that for any state which 
can be represented as a single valence bond state the en- 
tanglement entropy is equal to the number of EPR-pairs 
shared by the subsystems [30] (in our case this number 
is and 2 at nt s for the odd bipartitions and 1 for the 
even) . The oscillation follows closely, but not exactly, the 
curve 

<^)«2(l-cos 4 (g)). (9) 

The singlet and triplet populations at adjacent sites are 
quantities that can be experimentally probed via singlet- 
triplet spectroscopic measurements and Stern-Gcrlach 
techniques [11]. In terms of bond operators (see Eq. (4)) 
these quantities are defined as: 



j 

Seven(t) = ^ (f) |4j,2j+1 § 2j,2j+l i 

j 

W) - ^E^(*)I^A+ a ^+TA- +a ^(*)>. 

Sodd(t) = ^^2(lp(t)\slj + l,2j+2 S 2j + l,2j+2\^(t)) . 
j 

Using the evolution operator (6), the singlet-triplet 
populations can be shown to evolve as 

= \ ( 1 + 3cos4 (S)) ' (10) 

*£&'*(*) = s odd(t) = l - ■ 

The coherence of the singlet-triplet oscillations can 
help to characterize the quality of the dynamical evo- 
lution. These measurements, however, are only local and 
do not give any indication of the distance between the 
entangled atoms generated at t = (2n + l)t s . The latter, 
on the other hand, can be probed by measuring density- 
density correlations of the expanding cloud or noise cor- 
relations [16] 

eWM.W^I^^r-i, (ID 
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FIG. 3: The noise correlations for two periods in the sin- 
gle switch, t s = Tih/J. Numerical TEBD simulation for the 
infinite lattice. 



where n Q( r ) is the atom number operator for the compo- 
nent a at position r after time of flight. G(Q(r),Q'(r')) 
is directly related to the momentum-momentum correla- 
tions of the atomic cloud at the release time, t R . Deep in 
the Mott insulator regime G(Q(r), Q'(r')) can be rewrit- 
ten in terms of spin operators as 

G(q,t R ) = (^ tR )\^^e i9a(i ~ j) (l + S ^ S ^ t ^ 

= l -8 qfi +A{q,t R ), (12) 

where q = Q — Q' and a is the lattice spacing (Fig. 
1). While the first term in Eq. (12) reproduces the in- 
terference peaks at reciprocal lattice vectors character- 
istic of the Mott insulator state (due to the bunching 
of the bosons), the second term A(q,t R ) provides addi- 
tional information about the spin order in the system. 
For example, if the system is released exactly at times 
tfj i = nt s when it is in a valence bond state of length I 
(here I = 1,3), A(q,t R ) will exhibit spatial oscillations 
with periodicity dictated by the distance between the en- 
tangled atoms (see Fig. 3) 

A(q,t R = nt s ) = i[l + cos( ? aO]. (13) 

We note that the factor N in the denominator origi- 
nates from the short-range character of the interactions 
and therefore the entanglement is only shared between 
pairs. It limits the applicability of noise correlations as a 
suitable experimental probe in systems with large num- 
ber of atoms. However, the 1/iV factor should not be a 
problem in current II? systems with approximately 20 
atoms per tube [25]. 



V. PERIODIC SWITCH 

A. Generic case 

We now consider the iterative sequence of switching 
off and on the couplings J\ , J 2 every t s . One might 
think that if at time t s one reverses the couplings from 
J\ = 0, J2 = J to J\ = J, J2 = Q, the dynamics will just 
return the state into its original form, i.e. from Eq. (8) 
to Eq. (2); however, this does not happen. On the con- 
trary, as a result of the evolution under the swap oper- 
ators, atoms separated by four lattice sites become now 
connected by triplet valence bonds and so, at the time 
t = 2t s , the state evolves into a valence bond state with 
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FIG. 4: The noise correlations during the periodic switch, 
t s = nh/J. a) TEBD Simulation in the limit N -> 00. The 
fact that all the entangled pairs are of the same length is 
reflected in the periodic pattern, b) Exact Diagonalization 
on an open chain, 2iV = 14. The superposition of the triplet 
valence bonds with different lengths in the intermediate state 
around t/t s — N — 1 leads to a very weakly structured signal. 



6 



1 = 5 (Fig. 1), 

i^= 2t .) (2) =ni*w 5 ). (14) 

3 

The successive repetition of the switching procedure leads 
to the propagation of the entanglement across the lat- 
tice and after n switches, performed at times kt s (k = 
1, . . . n), one obtains entangled pairs with length 2n + 1. 

In the experimentally relevant case of an open chain, 
the sequential incrementation of the length of the en- 
tangled pairs is stopped when one member of the pair 
reaches the boundary of the lattice. The pair is then re- 
flected and continues moving through the lattice with its 
length remaining unchanged. Consequently, when after 
N — 1 switches the pair initially located at the center 
of the lattice reaches the boundary, a particular state 
that has the maximal possible length of entangled pairs 
is formed. While for an odd number of double wells it 
corresponds to a state with an EPR-pair connecting the 
edges of the chain, 

(N-l)/2 

\i>t=t a (N-l)) = l*Q,2JV-l)n l*2j-l,2JV-2j-l)l*2j,2.ZV-2j)> 
j = i 

for even N the maximal length of entangled pairs is I — 
2N-1, 

(JV-2) /2 

|^=t,(JV-l)) = II l*2j,2JV-2j-2)l*2j+l,2yY-2j-l) • 

Since the entanglement entropy of the state partitioned 
into its left and right half is simply given by the num- 
ber of EPR-pairs connecting the two parts [30], the 
state \ipt=t s (N-i)) has maximal entanglement entropy 
Sn = N. This growth of the entanglement for the case 
2N = 10 is depicted in Fig. 2. 

B. Implementation of remote entanglement 
protocol 

As we have seen, by applying the iterative swapping 
procedure to an open chain it is possible to engineer a 
state which has maximally separated entangled atoms 
and largest bipartite entanglement. Such a state can 
have relevant application in lattice-based quantum infor- 
mation proposals due to its large persistency of entangle- 
ment because in this case N qubits have to be measured 
to disentangle the state. The persistency of entanglement 
quantifies the robustness of the entanglement to noise. 
We remark that in this respect a cat state (macroscopic 
quantum superposition state e.g. TT • • - It) + I 4-1 

• • • II)) ) is fragile as even a single local measurement is 
sufficient to reduce it to a product state. The state we 
are engineering has persistency of entanglement as large 
as that of a cluster state, which is one of the key pre- 
requisites for using it as a one-way quantum computer 
resource [31]. 



Moreover, the \ipt=t s (N-i)) state is an eigenstate of the 
N-th switching operator, |V>t=(vV-i)0 = \ipt=Nt s ) and 
thus after 2N switches the state will be rolled back to 
the initial nearest-neighbor triplet-product state. This 
property can be useful for experimentally probing the 
state and quantifying the fidelity of the procedure. For 
example, by measuring the quality of the triplet product 
state after 2N switches one can get information about 
errors that occurred during the swapping process. 

We also note that even though we focused on the 
case of an initial array of triplet pairs, similar considera- 
tions hold if instead of triplets one starts with singlets or 
changes the sign of the coupling constants (as it would 
be in the case of fermions). 

In addition, our swapping procedure can also be used 
for transporting a particular state of an atom without di- 
rectly moving the particles: If we initially prepare all the 
atoms in the same state, say J., except for the atom at site 
i which we prepare in state t, after n periodic switches 
the state | will be transferred to the atom located at site 
i + n. 

The long-range entanglement produced by the switch- 
ing procedure can be experimentally probed by noise- 
correlation measurements. Although for finite lattices 
the expected ideal pattern of well-defined peaks at t — 
nt s (see Eq. (12)) changes to one with less regular struc- 
ture due to the distribution of different valence bond 
lengths, Fig. 4 shows that it still contains relevant in- 
formation such as the formation of well-defined peaks at 
q = and q = n/a when the distance between entangled 
atoms becomes maximal. 



C. Non-ideal conditions 

Up to this point we have assumed that Eq. (1) ac- 
curately describes the many-body dynamics. However, 
defects such as holes or doubly occupied sites will make 
this assumption invalid. 

We should emphasize that there is one particular con- 
dition which makes the entanglement generation possi- 
ble despite the presence of holes. Namely, this occurs 
when the single-particle tunneling time is engineered to 
be commensurate with t s . However, if this condition is 
not satisfied, in general the holes will hinder the gener- 
ation of long-distance entangled pairs and they should 
be suppressed for example by implementing additional 
filtering schemes such as the one proposed in Ref. [32]. 

Additionally, even though Eq. (1) was derived by tak- 
ing into account only virtual particle-hole excitations, 
real particle-hole excitations will certainly take place dur- 
ing the dynamical evolution. They would lead to oscilla- 
tions on top of the effective Hamiltonian dynamics with 
amplitude J/U and periodicity ~ h/U . Therefore, in or- 
der to efficiently average them out one has to work in 
the strongly correlated regime, i.e. with the condition 
Un, out *C U, though this implies smaller time scales for 
the dynamical evolution. In typical experiments, work- 
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ing in a parameter regime where particle-hole excitation 
effects arc negligible requires a superexchange coupling 
J I h of the order of 1 kHz (t s ~ 1 ms) and thus for a sys- 
tem with approximately 20 lattice sites, it will take about 
10 ms to generate entanglement between the atoms at the 
edges of the cloud. 

Another aspect of our procedure is that the long- 
distance entangled pairs are generated by switching the 
interactions at specific moments of time. In practice how- 
ever one always expects switching time uncertainties 5t 
and therefore the interval between consecutive steps will 
not be exactly t s but t s + St. Such inaccuracies will ac- 
cumulate and will degrade the quality of the final state 
exponentially with the number of lattice sites and the 
number of switches made during the process. Defining 
the fidelity of a state as T = \{ipt=nt a Wttnt s )\ 2 i where 
\ipt=nt ) an d \^t=nt s ) are the ideal and actual states gen- 
erated after n iterations, one can estimate the degrada- 
tion of fidelity using Eq. (6) 

_ r , n(St) 2 N, , . 

T ~ W - 1 I ) (15) 

where Tq = \(ip(t = 0)\U 3 l%2j+i)| 2 is the fidelity of 
the initial state. 



D. Entanglement purification 

To overcome all the limitations mentioned above one 
can combine our periodic switching scheme with entan- 
glement purification protocols. Starting from a large en- 
semble of generated low-fidelity entangled pairs, these 
protocols distill a smaller sub-ensemble which has suf- 
ficiently high fidelity. Entanglement purification can 
be implemented in a spin-dependent 2D superlattice 
as follows: after creating an array of ID independent 
chains along x-direction by suppressing tunneling along 
y-direction, one can use our procedure to generate many 
parallel long-distance entangled pairs within the ID 
chains, i.e. an atom at the site will be entangled 

with one at (i+l, j). Then tunneling along the x-direction 
should be inhibited and the following iterative procedures 
be applied: 

1. Lower the intra- well barriers along the ^/-direction 
of a spin-dependent superlattice, allowing only one 
of the species to tunnel [33]. This will introduce 
Ising-type interactions VJ ^'Sf^jSf^j+i between 
atoms at adjacent sites along y-axis and therefore 
will couple entangled pairs at (i, 2j) - (i+l, 2j) with 
pairs at (i, 2j + 1) - (i + 1, 2j + 1) respectively. 

2. Combine the Ising interaction with single-particle 
rotations, realized with the help of external mag- 
netic fields, to implement the C-Not gate required 
for the purification schemes described in Ref. [19]. 

3. Measure the spins at the (i,2j) and (i + l,2j) wells. 
If they turn out to be parallel, keep the correspond- 



ing pair at (i, 2j + 1) and (i + 1, 2j + 1), otherwise 
discard it. 

4. Release the measured atoms and merge the (i, 2j) 
and (i, 2j + 1) wells into a single one. Repetition of 
the above protocol will distill from the low-fidelity 
pairs the ones with higher fidelity. 

Let us now briefly discuss the experimental realizabil- 
ity of such purification protocols. To date, one of the 
main problems is the experimental implementation of 
step 3 due to the difficulty of measuring individual states 
at adjacent lattice sites. These atoms are separated by a 
distance of the order of an optical wavelength and there- 
fore diffraction fundamentally limits individual address- 
ability. One advantage of our scheme is that the atoms 
in a pair that should be measured are in general sepa- 
rated by many lattice sites, but nevertheless when the 
measurement is performed on one of the pairs, nearest 
neighbor atoms are still affected. One possibility to over- 
come this problem has been proposed recently in Ref. [34] 
where the use of nonlinear atomic response has been 
suggested for coherent optical far-field manipulation of 
quantum systems with resolution of up to a few nanome- 
ters. The implementation of the proposals of this kind 
in the controlled lattice environment may allow proof-of- 
principle experimental demonstration of quantum purifi- 
cation ideas. 



VI. HOMOGENEOUS SWITCH 

An interesting question which arises from the dynamics 
of the periodic switch is what happens with the quantum 
state if the double wells are no longer decoupled com- 
pletely, but instead there exists a finite tunneling between 
them. One expects that in this case the propagation of 
valence bond states will be suppressed after some period 
of evolution. To address this question, in this section we 
consider the case of a homogeneous switch J1/J2 = 1 
(case (3) in our classification), which formally can be 
considered as a particular case of quench dynamics: we 
prepare the system in a ground state of one Hamiltonian 
- a triplet (singlet) product state - and then suddenly 
change the quantum Hamiltonian to a new one - the 
isotropic ferromagnetic Heisenberg Hamiltonian -, which 
determines the subsequent evolution. 

In contrast to the periodic switch evolution, whose 
general characteristics are independent of the singlet or 
triplet nature of the starting state, the dynamics of the 
homogeneous switch is strongly affected by the symme- 
tries of the initial state. Consequently, we consider the 
cases when the initial state is in a triplet (singlet) config- 
uration separately. However, before starting the discus- 
sion we first provide a general overview of the dynamics 
of quantum quenched systems. 
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A. Quantum quench: general discussion 

The time evolution of a quantum state after a quan- 
tum quench has recently attracted a lot of theoretical 
[18, 35, 36] and experimental [25, 37, 38, 39] interest, 
in part due to the possibility of varying in real time 
the parameters of the optical lattice. For example, low- 
dimensional systems prepared in a gapless state and sub- 
sequently quenched into an insulator state have been ex- 
perimentally studied, addressing questions such as re- 
laxation to thermal states and collapse and revival ef- 
fects. The dynamics of exactly solvable models, e.g. 
an Ising chain [40, 41, 42], have also been the topic 
of investigation due to the fact that these systems sat- 
isfy many conservation laws which lead to non-trivial 
equilibration phenomena. Such behavior has been at- 
tempted to be explained in terms of a generalized Gibbs 
ensemble [43]. From the numerical side, recent advances 
in time-dependent density matrix renormalization group 
(DMRG) and TEBD methods [20, 21, 44] have allowed 
to study the quantum dynamics in bosonic and fcrmionic 
ID systems [45, 46, 47]. The numerical simulations seem 
to support the absence of thermalization, however, these 
methods are restricted to small and intermediate time 
scales. The case of the quench from the gapped phase into 
the critical regime has been studied using conformal field 
theory by P. Calabrese and J. Cardy [18, 35]. Numerical 
calculations [47, 48] support their results. The quench 
dynamics between gapped states can also be attacked 
using methods of exact solutions [49] and also demon- 
strate interesting dynamics associated with the absence 
of thermalization. On the other hand, more conventional 
approaches based on perturbative methods [50] and dia- 
grammatic expansions [51, 52] inevitably show dynamics 
associated with thermalization scenario. 

In the present work we adopt a numerical approach to 
deal with the quantum quench dynamics and postpone 
the analytical treatment for future publications. 



B. Initially prepared triplet state 

Let us first consider the case of the homogeneous 
switch dynamics when the initial state is a product of 
triplet states (Eq. (2)). In order to gain a general un- 
derstanding of this system, we note that while the ini- 
tial state has broken rotational and translational sym- 
metries, the Hamiltonian at t > (ferromagnetic Hciscn- 
berg) possesses both of these symmetries. Although its 
low-energy excitations are dominated by the spin-wave 
Goldstone modes corresponding to the broken continuous 
(rotational) symmetry, the quantum dynamics involves a 
bunch of highly excited modes which know nothing about 
the spontaneous breakdown of the continuous symmetry. 
We therefore face a dynamical competition between the 
initial state with broken symmetries associated with the 
initial condition on the one hand, and the whole spec- 
trum reflecting both of these symmetries on the other. 



As a result of this competition we expect the emergence 
of a complex magnetic state and the growth of the en- 
tanglement entropy. 

As we have pointed out, for the correct description of 
quantum dynamics it is not sufficient to rely on a low- 
energy effective theory because the details of the spec- 
trum can play a significant role. On the other hand, 
if we start with a state which involves many excited 
states, the characteristic features of the dispersion re- 
lation of the low-energy modes can be not so important. 
Also, quantities studied below are invariant under time- 
inversion symmetry and therefore the dynamics of our 
problem should have the same common features as that 
of the antiferromagnetic Heisenberg model. As a result, 
some common mechanisms should define the generic fea- 
tures of the quantum dynamics of these models. It has 
been pointed out recently [18, 35] that this generic be- 
havior can be understood in terms of classically moving 
quasiparticles [18, 36], whose transport correlations are 
bounded by the light cone {horizon effect). We interpret 
our results on the basis of these ideas. 



1. Entanglement 

We first focus on the evolution of the entanglement. 
The spatially anisotropic and weakly entangled initial 
state evolves into a highly entangled state with restored 
translational symmetry. This behavior is signaled by the 
growth of the entanglement entropy and the rapid de- 
cay of the oscillations between even and odd bipartitions. 
In Fig. 5 we plot the entanglement entropy of blocks of 
different sizes in a finite lattice. The plot shows that 
for short times, after the recovery of translational in- 
variance, the finite-block entanglement entropy exhibits 
linear growth. A saturation to a value close to the max- 
imal Si = I occurs for longer times. This is in agreement 
with results obtained with the use of conformal field the- 
ory [18] that predict a saturation value proportional to 
The growth of the entanglement limits the applicability 
of the numerical method (TEBD), as reasonable matrix 
dimensions (e.g. x — 1000) are only valid for weakly 
entangled systems (Si < log 2 (x) ~ 10). Consequently, 
it is impossible to verify the exact behavior of the en- 
tropy for large blocks. However, since in the intermedi- 
ate time regime the dynamically evolved state in finite 
lattices does not show significantly lower entanglement 
compared to an infinite system, to study this regime one 
can make simulations directly in the infinitely extended 
periodic system, where the translational symmetry can 
be exploited. This allows to reduce the computational 
cost by a factor of N compared to the finite-lattice sim- 
ulations [22]. 

We study the crossover that takes place from the 
'linear '-growth regime where Si — Soo (Soo stands for 
S^ en , S^ d , for even and odd I respectively), to a sat- 
uration towards a constant value. It is probed by the 
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FIG. 5: The entanglement entropy for the homogeneous 
switch, v a = Jn/2h. a) TEBD simulation for 2N = 36 with 
open boundary conditions. Si approaches the line I. b) The 
crossover from linear growth to saturation. Inset: deviation of 
the finite-block entanglement entropy from the infinite value. 
The crossover is well characterized by the saturation time de- 
fined by the spin-wave velocity, i* = — . 
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FIG. 6: The experimentally measurable triplet (t x ' v ' z , see 
text) and singlet (s) occupations at adjacent lattice sites, 
v 3 = Jn/2h. Exact diagonalization for 2N=22 sites. The 
equilibration of t z and t x ' y is not complete. The oscillations 
of even and odd bond correlations around the same value sig- 
nal the recovery of the translational symmetry. 



case is much slower than in this system. The reason for 
this is that in a one-dimensional lattice model the sharp 
crossover is smeared out by lattice effects (which explain 
why Si < Soo even for t < l/v s ) and, more importantly, 
by the non-linear dispersion relation. Due to the latter, 
particles moving slower than v s have to be taken into ac- 
count, what results in a slower saturation of Si(t) to a 
constant value at t > t*. 

While long-range effects at the 'horizon' are deter- 
mined by the 'fast' spin waves and the results from con- 
formal field theory are applicable, the slow quasi-particles 
will be of great importance for understanding the effects 
related to short-range phenomena. 



quantity (see Fig. 5) 

AS, = Soo- Si. (16) 

This crossover is a direct manifestation of the horizon 
effect. In the case of conformal invariancc, where rcla- 
tivistic dispersion relation ujk — v\k\ is assumed, the dis- 
tance between entangled atoms is always smaller than 
2vt. The entanglement grows linearly as long as the 
horizon is smaller than the size of the block. For the 
open chain considered here, with the block situated at 
one of the edges, the horizon has to be twice as large as 
the block length [48]. This allows to define a crossover 
time t* — l/v when Si[t > t*) becomes a constant [18]. 
Fig. 5 shows that using the spin-wave velocity of the 
Heisenberg ferromagnet, v s — Jir/2h, the crossover in- 
deed takes place around t* = l/v s . However, comparing 
results of Fig. 5 with the results of the quantum quench 
in the X JTZ-chain [48] , we find that the crossover in our 



2. Singlet-triplet population 

To study further the dynamical relaxation and the 
recovery of broken symmetries, we plot in Fig. 6 the 
singlet-triplet population at adjacent sites + The 
data are obtained by using an unbiased exact diagonal- 
ization technique (Lanczos algorithm [24]) on an open 
chain with 27V = 22 sites. After a certain time interval 
trelax one expects that the quenched initial values de- 
cay into a quasistationary regime, which is destroyed at 
t > t rec — AN/2v s due to quantum recurrence, when the 
edges of the expanding light cone begin to interfere [36] . 
From Fig. 6 we determine that v s t re i ax ss 5; the value 
v s trec = 22 corresponds to a conservative lower bound of 
the recurrence time. In the quasistationary regime only 
oscillations around an average value are observed. These 
oscillations can be associated with the finite bandwidth 
of quasiparticle energies [35]. We define average values 
for an arbitrary operator O in the quasistationary state 
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as follows: 

(0} gs = 



t rec trelax Jt 



dt(0(t)) 



(17) 



Comparing the x, y, ^-triplet populations, we see a ten- 
dency of relaxation towards proximate values: 



odd) is 
\ even/qs 
(t X - V ) 



0.334, 
0.343 , 
0.267, 
0.272. 



The relaxation of the singlet-triplet occupation numbers 
at even and odd bonds towards the same quasistationary 
values indicates the restoration of the translational sym- 
metry, also suggested by the entropy calculations. The 
difference between t z and t x,v of about 0.06 implies that 
the rotational symmetry is not completely restored. This 
difference is stable for various lattice sizes and choices of 
treiax (which by definition allows a certain freedom in its 
choice). This is a direct indication of missing thermaliza- 
tion in the quasistationary regime. The dynamic state 
does not fully reflect the symmetries of the Hamiltonian. 



3. Structure factor 

In Fig. 7 we plot the 'time-dependent' structure factor 
A(g,i), which is experimentally accessible by measuring 
the noise correlations (12). Alongside the persisting peak 
at q = 0, the picture shows the formation of a smooth 
peak at q ~ n/3a for all times t > t* , which signals the 
development of an unusual type of magnetic state. We 
checked that for small lattices (2N = 20) the peak is 
stable for v s t < 100. The height of the peak NA(q,t) is 
independent of the lattice size or the type of boundary 
condition, thus revealing the short-range nature of spatial 
correlations in the system. 

In order to explicitly study the relaxation of the cor- 
relation functions and to understand the origin of the 
incommensurate peak in the noise correlations, we plot 
in Fig. 8 the real-space correlation function 



G+-(i,t)= J2 Mst(t)sy(t)) 



\i-3\=l 



and the quantity 



Q + ~(l,t)= £ (-lYRc{S+(t)S-(t)), 



(18) 



(19) 



\i-l\=l 



which indicates that the translational symmetry is re- 
covered for long-range correlations. The most interesting 
effect we observe in the correlation functions is the sup- 
pression of the ferromagnetic (positive) nearest-neighbor 
correlations and the development of weak antiferromag- 
netic (negative) correlations for next-nearest neighbor 
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FIG. 7: The noise correlations for the infinite lattice, TEBD 
simulation, v s — Jn/2h. A broad peak at qa/n ~ 1/2 appears 
at v a t ~ 1, which shifts towards qa/n ~ 1/3 with the time 
evolution. 



sites. This is the origin of the incommensurate peak in 
the noise correlations (Fig. 7). The large-distance proper- 
ties of the correlation function do not contradict the pre- 
dictions of conformal field theory [35]. For instance, the 
correlations are ferromagnetic and change from their ini- 
tial values only when the 'horizon' of quasiparticle pairs 
l(t) = 2v s t passes, although in this case we find that the 
horizon is not absolutely sharp. It is important to notice 
that, although the horizon moves with constant speed, 
the intensity of the correlations decays fast with larger 
distances and the correlation length remains finite. 

The observed mixed correlations can be interpreted as 
an implication of energy conservation. At time t = the 
whole correlation energy is stored in the short-ranged 
triplets; at t > the action of the evolution operator 
leads to the formation of longer-distance singlets between 
spatially separated sites. This singlet component persists 
for longer times and leads to the appearance of the an- 
tiferromagnetic component in the spin-spin correlation 
function. Therefore the local redistribution of correla- 
tion energy, revealed in the partial AF correlation, is one 
possible explanation for the emergence of mixed correla- 
tions. 



C. Initially prepared singlet state 

In this section we study the case of the homoge- 
neous switch (with dynamical evolution determined by 
the Heisenberg chain, J\ = J2 = J > , Eq. (2)), but 
instead of starting from a triplet product state, we now 
start from a singlet product state, 



(20) 
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FIG. 8: Simulation of correlation functions for the infinite 
lattice using TEBD, v s = Jn/2h. a) Evolution of real-space 
correlation functions at fixed distances: the plot shows the 
tendency of longer-distance correlations to restore transla- 
tional invariance. b) Antiferromagnetic correlations at dis- 
tance 1 = 2, and rapidly decaying ferromagnetic correlations 
for larger distances. The straight line marks the horizon of 
quasiparticles moving with spin- wave velocity v s . The plot 
resolves magnitudes larger than 10 -7 . 



For bosonic systems this state can be experimentally re- 
alized by time evolution of the initial triplet product state 
in the presence of a magnetic field gradient [11]. This ini- 
tial state also corresponds to the decoupled double-well 
ground state of the respective fermionic system, though 
in this case the evolving Hamiltonian is the antiferromag- 
netic Heisenberg model instead of the ferromagnetic one. 
However, since the dynamical evolution is independent of 
the overall sign of the Hamiltonian, the results discussed 
in this section will also hold for the fermionic system. 

Unlike the case of initially prepared triplet state, here 
the spherical symmetry is not broken, and the popula- 
tions of the x, y and z components of the triplets are 
equal. From Fig. 9 we extract that 

(teven)qs = ^ odd )v ~ 0.125, 

(Seven)qs = (s dd)qs = 0.625. (21) 
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FIG. 9: The singlet and triplet populations for the initial 
singlet product state in the infinite lattice, TEBD simulation, 
v 3 = Jn/2h. The translational symmetry is recovered. 



These values are a direct consequence of the energy con- 
servation, 7jjy(H(t)) = 0.375. In Fig. 10 we study the 
spatial correlations. Fig. 10a) shows a rapidly developed 
broad antiferromagnetic peak in the noise correlations 
and weak incommensurate peaks at small wave vectors. 
These are due to large-distance spinon correlations, de- 
picted in Fig. 10b). The fact that the correlations remain 
negative after the spinon horizon passes (Fig. 10c) can be 
interpreted as a memory effect of the initial singlet state. 
Fig. 10c) shows, by investigating the quantity Q A (l,t) 
(see Eq. (19)), that the translational symmetry is recov- 
ered in the long-range correlation functions, as is the case 
also for the short-range singlet and triplet correlations. 

In general, the prepared singlet product state, due to 
its initial spherical symmetry, does not exhibit the strong 
mixing of anti- and ferromagnetic correlations, as the 
triplet state does. Although the observed spinon corre- 
lations are interesting from the theoretical point of view, 
their weak effect on the noise correlations is barely mea- 
surable experimentally We also note that the spinon 
correlations may disappear on large time scales which 
are inaccessible numerically. 



VII. CONCLUSION 

In this paper we proposed a novel protocol which cre- 
ates, from a system of two-component atoms initially pre- 
pared in an array of triplet (singlet) pairs on neighbor- 
ing sites, an array of long-distance triplet (singlet) pairs 
across the lattice. The method allows parallel generation 
of many entangled pairs, and can have relevant appli- 
cations for the implementation of quantum purification 
protocols in optical lattices. We also find that by ap- 
plying the iterative swapping procedure in an open chain 
one can engineer a state in which any atom located in the 
right half of the superlattice is entangled with an atom in 
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FIG. 10: Correlations for an infinite chain using TEBD simulations, v 3 — Jir/2h. a) The noise correlations for the system 
prepared in the singlet product state. Besides the strong antiferromagetic peak there are incommensurate branches for small 
q. b) The real-time correlation function at different moments of time. The correlations converge to an exponentially decaying 
antiferromagnetic behavior. For larger distances (I > 6), a staggered component centered around some finite negative value 
can be observed, c) The demonstration of how the longer-range correlations remain negative after the passing of the horizon 
(for even and odd distances) . 



the left half. This state has maximally separated entan- 
gled atoms and persistency of entanglement as large as 
that of a cluster state, which makes it suitable for being 
used as a component of a one-way quantum computer 
[31]. 

We also studied the evolution of an initial triplet (sin- 
glet) product state under a Heisenberg Hamiltonian. An- 
alyzing various observables we showed that while the 
long-range properties of the evolving state are in agree- 
ment with those predicted by conformal field theory, the 
non-universal short-range properties (e.g. the develop- 
ment of a magnetic state with mixed correlations), are 
not captured by such theoretical treatments [41] and have 
to be analyzed more carefully. They might be a manifes- 
tation of a special type of thermalization (in the sense of 
generalized Gibbs ensemble [43]), observed in integrable 
systems. 

The analysis presented in this paper demonstrates that 
the coherent evolution of an initial state, which itself can 
be easily prepared - in our case it is just an array of 



triplet (singlet) states on neighboring sites -, is a feasi- 
ble way to generate complex magnetic states with cold 
atoms. The dynamical generation method is not con- 
strained by the difficulty of actual (physical) engineering 
of exotic Hamiltonians or by the low temperatures re- 
quired to reach their ground states. On the other hand, 
without a careful analysis it is difficult to predict a priori 
the properties of the non-equilibrium state into which the 
system evolves as a result of coherent quantum dynamics. 
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